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Abstract. We have investigated whether the Vi feature of NHj is a viable candidate for the 6.85 ^m absorption 
band seen towards embedded young stellar objects. To produce NH^ astrophysical ice analogs consisting of H2O, 
CO2, NH3 and O2 were UV photolysed. The IR spectra reveal peaks that are identihed with the NH^, NO^, NO;^ 
and HCO;^ ions. It is shown that the NHj matches two absorption features that are observed towards embedded 
young stellar objects, i.e., the strong 6.85 fim feature and the 3.26 fim feature. The characteristic redshift with 
temperature of the interstellar 6.85 fim feature is well reproduced. The abundance of NHj in interstellar ices would 
be typically 10 % relative to H2O. The experiments show that the counterions produce little distinct spectral 
signature but rather a pseudo-continuum if a variety of them is present in a H2O dominated environment. The 
anions could therefore go undetected in IR spectra of interstellar ice. In the ISM, where additional mechanisms 
such as surface chemistry and additional elements such as sulfur are available many acids and an even wider 
variety of anions could be produced. These components may be detectable once the ices sublime, e.g, in hot cores. 

Key words. Methods; laboratory - Stars: individual: W33A - Stars: individual: MonR2:IRS3 -ISM: abundances 
- ISM: molecules - Infrared: ISM: lines and bands 



1. Introduction 

The nature of the 6.85 /im absorption feature towards 
embedded Young Stellar Objects (YSO's) has remained 
an enigma since its discovery 25 years ago (Russell 
et al. 1977). While the first high resolution observa- 
tions of this band were recently obtained by the Short 
Wavelength Spectrometer (SWS) on board the Infrared 
Space Observatory (ISO; Schutte et al. 1996, Dartois et 
al. 1999a; Keane et al. 2001), this did not yet help to 
clarify its origin. A number of candidates, such as carbon- 
ates, and the CH deformation modes in organic molecules 
like methanol, could be excluded based on the feature's 
spectral properties and the absence of strong additional 
bands. One possibility, however, the z/4 mode of the am- 
monium ion (NH|), first proposed by Grim et al. (1989b), 
could not be excluded. NH4 is produced in astrophysi- 
cal ice analogs by acid-base reactions. This is achieved 
either by deposition and warm-up of NH3 together with 
acids such as HNCO or HCOOH (Novozamsky et al. 2001, 
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Schutte et al. 1999), or by photolysis of ice mixtures con- 
taining NH3 (e.g., NH3/O2; Grim et al. 1989a) where the 
acids are produced by the photochemistry. Recently, the 
good spectral correspondence between the 1^4 feature of 
NH^ and the interstellar 6.85 /xm absorption as seen by 
ISO/SWS was again demonstrated (Demyk et al. 1998). 
Nevertheless, this assignment faces a fundamental prob- 
lem, since it requires a large abundance of NH^ of ~ 10 
% in the ices near YSO's. An equal amount of negative 
charge would need to be present. Although the ions OCN~ 
(also referred to as NCO- in the chemistry literature) and 
(probably) HCOO~ have been identified, their abundance 
falls far short of what is needed for the balance (Schutte 
et al. 1996a; Gibb et al. 2000; Keane et al. 2001; 2002). 
No sign of other negative solid state species has shown up 
in the ISO data. 

This paper revisits the possibility that NH^ is respon- 
sible for the interstellar 6.85 /im band. We investigate in 
the laboratory whether a variety of anions can supply the 
required countercharge, and that such a combination of 
species is either lacking in outstanding IR features, or fall 
below current detection limits. Taking the cosmic abun- 
dances of the elements into consideration, the counterions 
would have to consist of H, O, C, N and S in order to make 



W. A. Schutte et al.: NH4 and the 6.85 fim band near Young Stellar Objects 



2 

a significant contribution. If we furthermore constrain the 
acid precursors to species with at most 1 carbon atom, 
possibilities would still include a variety of thermodynam- 
ically stable species, i.e., OH", CN", OCN", HCOQ-, 
HCO3 , CO^", NO", NO", SO^-, HSO3 , and SO^^. In 
principle an astrophysical ice analog containing NH^ and 
a mixture of such negative ions could be prepared by de- 
positing NH3 together with a variety of acids. However, 
experimental limitations compel us to produce the ions 
in situ with UV photolysis. To this end, we processed 
mixtures of H2O/CO2/ NH3/O2 to produce acids such as 
H2CO3 (carbonic acid), HNO3 (nitric acid) and HNO2 (ni- 
trous acid). These acids react with NH3 forming NH4 and 
a mixture of counterions, as desired. Earlier work on sim- 
ilar mixtures indeed showed the formation of such species 
(Moore & Khanna 1991; Gerakines et al. 2000; Grim et 
al. 1989b). Of course, this experimental method of pro- 
ducing the ions may be different from what happens in 
the ISM. For example, the high abundance of O2 in our 
samples, necessary to stimulate the photochemical pro- 
duction of oxygen rich acids, may considerably exceed the 
abundance in interstellar ices (Vandenbussche et al. 1999). 
In space acids are possibly produced by other mechanisms 
such as surface chemistry (e.g., Keane & Tielens 2002, in 
preparation). 

To verify whether NH^ is a plausible candidate for 
the identification of the 6.85 /im band a number of cru- 
cial issues need to be adressed. First of all, detailed spec- 
troscopy of the 1^4 mode is needed to investigate its be- 
haviour as a function of ice composition and tempera- 
ture. Next, the issue whether other infrared features of 
NH^" could be present in interstellar ice spectra must be 
adressed. Subsequently, the counterions that are produced 
in the experiments should be identified to assess their as- 
trophysical relevance. Finally, the spectral signature of the 
counterions under various conditions is studied, where the 
prime question is whether conditions exist at which their 
infrared signatures become inconspicuous. 

The paper is organized as follows. In Sect. 2 we review 
the experimental techniques. Sect. 3 summarizes the re- 
sults. NII4" is created by UV photolysis of astrophysical ice 
analogs. We describe the spectral properties of NHJ under 
such conditions. Furthermore we establish which counte- 
rions are formed in the experiments and study their spec- 
troscopic properties as well. In sect. 4 the obtained NH^ 
spectra are compared with observations of YSO's. Besides 
matching the 6.85 fim band, it is furthermore shown that 
the interstellar 3.26 /im feature closely corresponds to one 
of the NH^ features. Also, using the experimental results, 
it is investigated whether the absence of features due to 
counterions in the observations can be reconciled with an 
assignment of the 6.85 /xm band to NH^. In Sect. 5, the 
astrophysical implications of the NH^ identification are 
discussed. Sect. 6, finally, summarizes the conclusions of 
this paper. 



2. Experimental 

Detailed descriptions of the general procedure for creation 
and photolysis of ice samples and the measurement of their 
infrared spectra have been published earlier (Gerakines et 
al. 1995; 1996). In summary, the set-up consists of a high 
vacuum chamber (10~^ mbar), with an IR transparent Csl 
substrate mounted on a cold finger which is cooled to ~ 
12 K. Samples were slowly deposited (~ 3 x 10^^ molec. 
cm^^ s^^/4 /zm hr~^) through a narrow tube controlled 
by a regulation valve. Photolysis by vacuum UV was sub- 
sequently performed by a hydrogen flow discharge lamp 
('^ 5 X lO^'' photons cm^^ s^^; F,photon > 6 eV). For 
thorough photolysis, the thickness of the sample should 
be '~< 0.2 fim. To overcome this limit, in one case depo- 
sition and photolysis were performed simultaneously. In 
this way a photolysed sample can be produced of several 
microns thickness, which greatly enhances the amount of 
photoproducts and thus the S/N ratio of the IR spectrum. 
Afterwards the sample was warmed in steps. The evolu- 
tion of the sample throughout the photolysis and warm-up 
sequence was monitored by infrared transmission spec- 
troscopy. In addition to the photochemical experiments, 
we did band strength measurements for NH^ by deposi- 
tion and warm-up of ices containing NH3 and HCOOH. 

The reagents used in these experiments were II2O (pu- 
rified by three freeze-thaw cycles), CO2 (Praxair, 99.996% 
purity), NH3 (Praxair, 99.99 % purity), O2 (Praxair, 
99.999 % purity) and HCOOH (Baker, 98 %). When 
ice samples containing both NH3 and CO2 or NH3 and 
HCOOH were prepared, the CO2 or HCOOH was de- 
posited through a separate tube, to prevent reactions with 
NH3 prior to deposition. This method is described in 
Gerakines et al. (1995). 

Table 1^ provides a log of our experiments. Column den- 
sities and abundances of H2O, CO2 and NH3 were directly 
obtained from the IR spectrum using band strengths from 
the literature (Gerakines et al. 1995; Kerkhof et al. 1999). 

3. Results 

3.1. Photolysis of R2O/CO2/NR3/O2 mixtures 

We photolysed a number of mixtures of 
H2O/CO2/NH3/O2 to produce NH| and a variety 
of negative ions (Table |l|). Fig. [l] shows the evolution of 
spectrum of photolysed H2O/CO2/NH3/O2 = 10/2/1.1/1 
(Expt. 4; Table |l|) during warm-up from 12 K to 220 K. 
First of all we want to verify whether NH^ is formed. 
Therefore, Fig. ^ also shows the spectrum of the binary 
mixture HNCO/NH3 = 1/1.2, after warm-up to 120 K 
(no photolysis). In this sample NH^ is readily formed by 
proton transfer between the isocyanic acid and the am- 
monia (Novozamsky et al. 2001). The NH^ features fall 
at 3000 cm-i (broad; vi), 3200 cm^^ {21^2), 3060 cm^^ 
{lyi+iys), near 2860 cm^^ {^i^a) and near 1450 cm^^ (^^4)1 
where is a lattice mode (assignments from Nakamoto 
1972). The first four features blend together in a broad 
structure which extends from ~ 3500 - 2400 cm~^. 
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Table 1. Log of experiments 







Mixture 




N 


d" 


dose 


dose 




H2O 


CO2 


NH3 


O2 


molec. cm^^ 


fim 


photons cm~^ 


eV molec.-" 


1 


10 


10 





10 


7.2(17) 


0.37 


9(17) 


<11= 


2 


10 


5 


7 


10 


4.0(17) 


0.17 


9(17) 


<20= 


3'' 


10 


2.2 


1.2 


1 


2.5(19) 


9.5 


2.4(19) 


8.6 


4 


10 


2 


1.1 


1 


8.5(17) 


0.32 


6(17) 


<6.4" 


5 


10 


1.2 


0.9 


0.9 


4.8(17) 


0.17 


9(16) 


<1.7" 


6 





1.2 


1 


1 


4.4(17) 


0.23 


9(16) 


<1.8'= 


7 








1 


1 


4.5(17) 


0.18 


9(17) 


<18'^ 



Thickness obtained from the column density using a nominal density of 1 g cm 
Obtained by simultaneous deposition and photolysis. 
Sample not optically thick, therefore dose is an upper limit. 
Using an average photon energy of 9 eV (Jenniskens et al. 1993). 
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Fig. 1. Spectra of the photolysed sample 
H20/C02/NH3/02= 10/2/1.1/1 (expt. 4): a. Before 
warm-up (12 K); b. at 120 K; c. at 220 K. Curve c was 
scaled up by a factor 3. Curve d shows the spectrum of 
HNCO/NH3 = 1/1.2 at 120 K (after the transformation 
of most of the original species into OCN" and NH^; 
arrows identify features of NH^). The spectra have been 
offset for clarity. 



Directly after the photolysis of the H2O/CO2/NH3/O2 
mixture, the feature is already clearly present at ^ 
1500 cm-^ The H2O ice absorption in the 3400- 2800 
cm"^ region obscures the other absorptions due to 
Upon warming to 220 K, causing the ice to sublime, 
the other NH4" features can all be seen. Clearly, the 
ammonium ion is produced by the photolysis. 

As discussed in Sect. 1, we want to investigate the spec- 
tral properties of the z/4 NH^ feature. To study the tem- 
perature dependence. Fig. |2| shows the 1900 - 1000 cm"^ 
spectra of photolysed H2O/CO2/NH3/O2 = 10/2/1.1/1 
(expt. 4). Spectra are shown directly after photolysis at 
12 K, and at a variety of temperatures up to 240 K. Some 
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Fig. 2. Thermal evolution of the photolysed ice mixture 
H20/C02/NH3/02= 10/2/1.1/1 (expt. 4), blow-up of the 
1900 - 1000 cm^^ region: a. before warm-up (12 K); b. 
after warm-up to 120 K; c. to 150 K; d. to 180 K, e. to 
220 K; f. to 240 K. Apart from the multiplication factors 
indicated in the figure, the vertical scale is identical for all 
curves. However, the spectra have been offset for clarity. 



important spectral properties can be gleaned from the fig- 
ure. Upon warm-up to 120 K the depth of the feature 
near 1460 cm^^ increases by a factor ~ 2. While this in- 
crease could result from further acid-base reactions dur- 
ing the warm-up, the absence of any feature attributable 
to acids such as H2CO3 or HNO3 after the photolysis 
(Fig. |, H2CO3 has bands at 1727, 1480 and 1275 cm"^; 
Gerakines et al. 2000; HNO3 has a strong feature at 1300 
cm"^; McGraw et al. 1965) argues against this possibil- 
ity. Thus, the growth of the NH^ feature during warm-up 
is likely caused by an increase of its intrinsic strength. 
Furthermore, the feature shifts strongly redward with 
temperature. The enhancement and the shift reflect the 
strong interaction of the ion with its environment. Fig. |^ 
plots the position of the i^^ band vs. temperature for a 
number of samples (expts. 2, 4 and 6 of Table |l|). It can 
be seen that, while there is a steady redshift throughout 
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Fig. 3. Redshift of the i'4NH^ feature as a function of tem- 
perature: H2O/CO2/NH3/O2 = 10/2/1.1/1 (expt. 4; open 
circles); 1/0.5/0.7/1 (expt. 2; open squares); 0/1.2/1/1 
(expt. 6; filled circles). 
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Fig. 4. The thermal evolution of the residue of photol- 
ysed NH3/O2 = 1/1 (Expt. 7). From central to top: the 
residue at 180 K; at 220 K; subtraction of the 220 K from 
the 180 K spectrum. The vertical scale is identical for 
all curves, however, the spectra have been offset for clar- 
ity. For comparison, ammonium nitrate (NH4NO3) and 
sodium nitrite (NaN02) are shown at the bottom (Miller 
& Wilkins 1952). 



the warm-up, the shift is particularly pronounced during 
H2O sublimation between 160 and 180 K. 

As discussed in Sect. 1, it is essential to identify the 
negative ions produced in the experiments to assess their 
astrophysical significance. To this end, we study the IR 
spectrum of the residue which remains after the ices have 
fully evaporated. As a basis for the interpretation of more 
complex samples, Fig. ^ shows the warm-up behaviour 
of the residue of photolysed NH3/O2 = 1/1 (expt. 7; 
Table |l|). The evolution is best understood by a "hot to 



cold" interpretation, i.e., starting with the most refrac- 
~ tory component at the final stages of the warm-up, and 
then working backwards in the sublimation sequence to 
see which spectral elements are sequentially lost. After 
warm-up to 220 K, a component remains which is stable 
up to ^ 280 K. Its spectrum resembles that of ammonium 
nitrate (Fig. ^; NH4NO3; throughout the paper we will 
use the standard chemical notation for salts. It must be 
noted however that the constituents of these compounds 
_ are ions, in this case, NHJ and NO3"). The NO;^ ion was 
previously identified in photolysed NH3/O2 ices (Grim et 
2 al. 1989b). However, while there is a general correspon- 
dence with the literature spectrum of ammonium nitrate 
there are some differences in the relative intensities and 
positions of the bands. The literature spectrum was pro- 
duced from fine crystalline powder in nujol mull, while our 
samples consist of a mixture of salts at low temperature. 
Spectra of ions are generally quite sensitive to factors like 
temperature, matrix and degree of annealing, due to the 
strong interaction of the ions with the matrix (c.f. OCN^ 
in various salt matrices; Maki & Decius 1959). This ef- 
fect possibly causes the differences. During warm-up from 
180 to 220 K two strong features disappear at 1270 and 
1220 cm~^, together with a fraction of the 1450 cm^^, and 
2400 - 3500 cm-i complex due to NH^. The band at 1220 
cm^^ may correspond to the main feature of NOj^. As 
compared with the literature spectrum of NaN02 (in nu- 
jol mull; Fig. ^), the band is shifted by 35 cm~^, possibly 
due to matrix interactions. The assignment of this com- 
ponent with ammonium nitrite (NH4NO2) is supported 
by the identification of the NO2 ion in the photolysed ice 
(Grim et al. 1989b). The 1270 cm^^ feature corresponds 
to an unidentified photoproduct. 

Fig. H shows the spectral evolution of the residue of 
photolysed H2O/CO2/NH3/O2 = 10/2.2/1.2/1 (expt. 3). 
To probe the nature of the anions, we will again analyse 
this diagram "hot to cold". As with the NII3/O2 sample, 
the most refractory component, remaining at 280 K af- 
ter all other material has sublimed, resembles ammonium 
nitrate (Fig. ||). During warm-up from 240 to 280 K a 
component sublimes which is characterized by two strong 
bands at 1600 and 1300 cm~^, together with four weaker 
features at 1375, 1010, 830, and 690 cm^^. The disap- 
pearance of these bands is accompanied by a decrease of 
the NH^ features. All these features have counterparts in 
the spectrum of ammonium bicarbonate in a KBr pellet 
(NH4HGO3; Fig. I). StiU, some differences are apparent in 
width and relative intensity, while furthermore some extra 
bands are present in the literature spectrum of which only 
minor indications are seen in the 240 to 280 K component. 
Again, we note that the spectra of ions are generally quite 
sensitive to such factors as temperature, matrix and de- 
gree of annealing, and it seems likely that the differences 
could derive from such factors. 

To verify the presence of HCO3" in the 240 - 280 K 
fraction, we photolysed the ice mixture H2O/CO2/O2 = 
1/1/1 (Expt. 1, Table |l|) to see whether its precursor, 
carbonic acid (H2CO3) is formed when NH3 is omitted 
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Fig. 5. The thermal evolution of the residue of photolysed H2O/CO2/NH3/O2 = 10/2.2/1.2/1 (Expt. 3). From bottom 
to top: The spectrum of the sample at 200 K; at 240 K; at 280 K; subtraction of the 240 K from the 200 K spectrum; 
subtraction of the the 280 K from the 240 K spectrum. The vertical scale is identical for all curves, however, the spectra 
have been offset for clarity. For comparison, ammonium nitrate (NH4NO3; Miller & Wilkins 1952) and ammonium 
bicarbonate (NH4HCO3; Khanna & Moore 1999) are shown. 



from the mixture. The photolysis gave rise to features at 
3040, 2840, 2627, 1727, 1480, and 1275. These bands can 
all be ascribed to carbonic acid (Moore & Khanna 1991; 
Gerakines et al. 2000). 

The residue fraction subliming between 200 - 240 K is 
characterized by a rather complex spectrum. Features at 
1605, 1375, 1300, and 1010 cm"! show that some HCO3 
sublimes in this temperature interval. However, the weak- 
ness of these bands shows that the contribution of bicar- 
bonate to the total negative charge in the 200 - 240 K 
fraction is minor. The band at 1670 cm~^ and the broad 
structure between 1000 and 600 cm~^ likely correspond 
to a small quantity of H2O which is embedded in the 
more refractory material (Schutte & Buys 1961; Ryskin 
1974). Likewise, the 2340 cm~^ feature can be ascribed 
to a small amount of embedded CO2. Like for the pho- 
tolysed NH3/O2, the 1220 cm^i may be due to the NO^ 
ion (Fig. However, the intensity of this band shows 
that NO^ accounts for only ^ half of the countercharge 
of the NH^ which sublimes in the 200 - 240 K inter- 
val. The 2170 cm~^ feature indicates the production of 
some OCN~ (Hudson et al. 2001), however, it is too weak 



to make a substantial contribution to the countercharge. 
We conclude that an important fraction of this relatively 
volatile residue component must consist of unknown an- 
ions. No sign of such species, or of the acids from which 
they originate, was found in experiments with less than 
four initial ice components (i.e., H2O/CO2/O2, NH3/O2; 
sect. 3.1). It therefore seems probable that they are rela- 
tively complex ions whose formation pathway involves all 
four original components. 

As outlined in Sect. 1, a prime issue for the viabil- 
ity of the NH^ assignment of the 6.85 /im band concerns 
whether it is possible to produce a clear spectral sig- 
nature while features of counterions are weak or absent 
as for the interstellar spectra. From Fig. ^ it can be seen 
that intense broad bands near 1320 and 1220 cm~^ are 
present once the ices have sublimed at ^ 180 K. As dis- 
cussed above, these bands are caused by negative ions, 
probably HCO3 , NO;^, and NO^. Before ice sublimation, 
these bands are however quite inconspicuous, resulting in 
a smooth sloping spectnun in the 1400 - 1200 cm~^ range 
with only minor substructure due to the 1300 and 1270 
cm~^ features of unidentified photoproducts. 
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Fig. 6. Spectra of photolysed ice samples, after warm- 
up to 160 K, enlargement of the 1900 - 1000 cm^^ re- 
gion: a. NH3/O2 = 1/1 (expt. 7); b. HsO/COa/NHs/Os^ 
1/0.5/0.7/1 (expt. 2); c. 10/2/1/1.1 (expt. 4); d. and 
10/1.2/0.9/0.9 (expt. 5). The spectra have been offset for 
clarity. 



To further investigate the influence of the H2O ma- 
trix on the appearance of the anion bands, Fig. ^ com- 
pares the 1900 - 1000 cm~^ spectra of the ice samples 
H2O/CO2/NH3/O2 0/0/1/1, 1/0.5/0.7/1, 10/2/1/1.1 
and 10/1.2/0.9/0.9 (Expts. 2, 7, 4 and 5) AU spectra were 
taken after warm-up to 160 K and sublimation of NH3, 
which enables a better view of the features of the pho- 
toproducts. Samples 2 and 7 show strong bands due to 
negative ions at 1385 and 1335 cm"^ (NO;^) and at 1230 
cm^^ (NO2'). There is hardly a sign of these or other fea- 
tures of anions in the H2O dominated samples (expts. 4 
and 5), even though the intensity of the NH^ band at 
~ 1475 cm~^ is similar. This again shows that in the H2O 
dominated samples features due to anions are inconspicu- 
ous. 

3.2. Summary 

The photolysis of H2O/CO2/NH3/O2 ice mixtures effi- 
ciently produces NH^ . These experiments revealed the fol- 
lowing four important spectral properties of NH^ and the 
negative counterions which pose important constraints for 
an identification of NHj in interstellar ices: 

1. During warm-up from 12 to 280 K, the NH4" 1/4 
feature shows a pronounced shift from ^-^1480 cm^^ to 
-1420 cm-i. 

2. Between 12 to 120 K the NH^ 1/4 feature grows 
smoothly by a factor 2. 

3. In H2O ice three relatively weak features of NH^ are 
only evident after the ice sublimation at ^ 180 K. These 
are the 2i/2 band near 3200 cm~^, the vi+v^ band near 
3060 cm~^, and the 2v4 band near 2860 cm~^. 

4. In H2O dominated ice, the spectral features of the 
counterions are only apparent after the sublimation of the 



Table 2. Band strength of the V4 NHj band 



mixture 


T (K) 


A (cm ion ^) 


H2O/NH3/HCOOH 






100/3.6/3.6 


120 K 


4.4 (-17) 


NH3/HCOOH=4/10 


80 K 


4.0 (-17) 


H2O/CO2/NH3/O2 






10/1.2/0.9/0.9, + hv 


150 K 


>2.7 (-17) 



ice, but are unconspicuous in the II2O dominated ice ma- 
trix. 

Concerning this last point, analysis of the spectrum of 
the residue shows that HCO;^, NO;^, NOj^ are probably 
present. The total abundance of these species is somewhat 
unsufficient for balancing the NH^, indicating the produc- 
tion of additional, yet unidentified, anions. 

The identifications of anions from the IR spectra are in 
some cases tentative, while in other cases no identification 
was possible at all. A positive identification of the anions 
and corresponding acid neutrals may help to constrain 
the composition of interstellar ices (see Sect. 5 below). In 
future, the nature of the negative species could be studied 
by analysis of the molecules that evaporate during the 
warm-up, for example, by mass spectroscopy. 

3.3. Strength of the U4 NH^ band 

To obtain column densities, it is necessary to measure the 
intrinsic strength of the NH^ 1/4 feature. This can be de- 
rived from experiments in which NH^ is formed by simple 
warm-up (no photolysis) of an ice containing NII3 and 
an acid. In this case the quantity of NH^ formed dur- 
ing warm-up is equal to the amount of NH3 which is lost. 
Since the band strength of the NH3 umbrella mode in var- 
ious matrices is known (Kerkhof et al. 1999), this allows a 
precise determination of the NII4" band strength. We per- 
formed two experiments involving NII3 with formic acid 
(HCOOH; see Schutte et al. 1999 for experimental details). 
The results are listed in Table Since the strength may 
depend on the composition of the matrix and nature of the 
counterion, the band strength of NHj in the photolysis 
experiments may differ from the IICOOII/NII3 mixtures. 
No direct measurement in the photolysis experiments is 
possible, because, besides acid-base reactions, other pro- 
cesses may contribute to the NH3 destruction during the 
photolysis. Therefore the amount of NH3 destroyed ex- 
ceeds the formation of NII4" by an unknown factor and 
only a lower limit to the NHJ band strength can be ob- 
tained. The most constraining lower limits are derived for 
the experiments with H2 0-dominated ices (see Table |). 
This is not surprising, since in these experiments features 
of N-containing photoproducts other than NHj are small 
(Sect. 3.1). 

We will adopt a standard band strength A(i/4,NII^) 
— 4.4x10^^^ cm molec.^^. This is close to the value of 
(2.5 - 3.5) X 10^^^ cm molec.^^ found in aqueous solution 
(Lowenthal & Khanna, in preparation). As discussed in 
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Sect. 3.1, the band strength varies by a factor 2 with 
temperature in our photolysis experiments. It is unclear 
whether the standard band strength corresponds to the 
highest or lowest bandstrength in the photolysis experi- 
ments. We will therefore adopt a standard uncertainty of 
a factor 2 in this value. 



4. Comparison to observations of Young Stellar 
Objects 

In this section we investigate a number of criteria for the 
identification of the interstellar 6.85 /im absorption with 
the V4 feature of the ammonium (NH^) ion. First of all, 
the feature should provide a good match to the observed 
band. Second, the characteristic redshift of the feature 
with temperature (Sect. 3.1; Fig. |^) should show up in 
the observational data. Third, additional bands of Nllt 



X (|im) 



are sought, specifically the 2^2, vi + vz, and 21^4 features 
near 3200, 3060, and 2860 cm^^ (3.12, 3.27, and 3.50 
/im). Fourth, there should be no spectral structure due 
to counter-ions that is inconsistent with the interstellar 
data. It will be shown that the available observational in- 
formation indicates that all these criteria are satisfied. For 
further reference. Table ^ lists the relevant observational 
data for all lines of sight where high quality spectra of the 
6.85 /im feature (i.e., by ISO-SWS) are available. Data on 
the interstellar 3.26 and 3.48 /xm features are included be- 
cause they are possibly associated with the NH4 vi -\- 
and 2i/4 bands (see below). The depth of the 6.85 /im 
band was measured from the original data. These were ob- 
tained from the ISO-SWS database (for NGC7538:IRS1) 
and from Keane et al. 2001 (all other objects). The depth 
was obtained by subtracting a linear baseline in the log(F) 
vs /im plane. The baseline was drawn through the 5.5 and 
7.5 /im points. 

4.1. Comparison with the interstellar 6.85 fim feature 
4.1.1. W33A 

Fig. I shows the 4000 - 1200 cm'^ (2.5 - 8.3 /im) spectra 
of the high mass YSO W33A obtained by ISO/SWS in ob- 
servation modes 1 and 6 (de Graauw et al. 1996). Spectra 
were adopted from Gibb et al. (2000) and Keane et al. 
(2001), respectively. While the SWS6 spectrum gives the 
highest S/N, the wide range of the SWSl spectrum allows 
an accurate baseline determination. We defined the base- 
line over the entire 4000 - 1200 cm~^ range by drawing 
a smooth third order polynomial through the continuum 
regions 4100 - 3700 cm^^, 2010 - 1980 cm~\ and 1840 - 
1810 cm~^ (Fig. 0). This procedure differs in some respects 
from the baseline definition in Keane et al. and Gibb et al. 
Keane et al. defined the baseline only over a limited data 
range (1900 - 1300 cm~^), and included the region around 
1300 cm^^, just before the onset of the silicate feature, in 
the baseline fit. It is however clear from Fig. |^ that, at 
least for W33A, there is considerable absorption in this 
region. Gibb et al., while also fitting the baseline over a 



E 




Fig. 7. ISO/SWS 1 (thin line) and SWS 6 spectra (1900 
- 1200 cm-i thick fine) of W33A (from Gibb et al. 2000 
and Keane et al. 2001, respectively). The smooth dashed 
curve is a third-order polynomial continuum fit. 
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Fig. 8. Correction of the 6.85 /im band of W33A for the 
contribution of CH3OH. a. H2O/CH3OH/CO2 = 1/1/1, 
at 112 K (from Ehrenfreund et al. 1999; band strength has 
been scaled to correspond to the CII3OH column density, 
see text); b. Optical depth spectrum of W33A obtained by 
subtraction of the continuum from the SWS 6 data (See 
Fig. 0); c. Corrected optical depth spectrum of W33A ob- 
tained by subtraction of curve a from curve b. The vertical 
scale is identical for all curves, however, curve c has been 
offset for clarity. 



similar broad range, included the 2500 - 2440 cm^^ region 
in the fitting procedure. However, it is probable that both 
the long- wavelength shoulder of the II2 O band ( Willner et 
al. 1982; Smith et al. 1989, see also Fig. |^) as well as the 
broad II2O combination feature near 2220 cm~^ (Gibb et 
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Table 3. Properties of the 3.26, 3.48, and 6.85 /im features towards YSO's. 



Object 


N(H20)'' 


F(45) b 
F{100) 




3.26 fim" 


3.48 fim" 


6.85 




t(3.26) ^ 
t(6.85) 


IN t^iN tl^ } 

N(H20) 




lO^^cm"^ 


u{cm ^ 


) 


T 






% 


W33A 


12 


1.27 




<0.4 


>0.29 


1471 


1.1 


<0.33 


11 


NGC7538:IRS9 


7 


2.07 






0.13 (0.02) 


1473 


0.29 




5 


GL2136 


5.0 


2.75 




<0.14 


0.14 (0.02) 


1459 


0.23 


< 0.6 


6 


S140:IRS1 


2.8 


3.24 


3072 


0.036 (0.007) 


0.027 (0.007) 


1451 


0.16 


0.23 (0.06) 


7 








3085S 


0.050f 


0.050" 






0.31 




W3:IRS5 


5.8 


4.12 






0.13 


1459 


0.25 




5 


MonR2:IRS3 


1.6^ 


4.32 


3071 


0.049 (0.007) 


0.036 


1437 


0.23 


0.21 (0.03) 


17 


AFGL7009S 


12 










1462 


1.1 




10 


GL989 


3 




3098'' 


0.069» 


0.069" 


1471 


0.11 


0.63 


4 


p Oph Elias 29 


3 






<0.03 


0.089 (0.007) 


1470 


0.07 


<0.4 


3 








3080S 


0.085» 


0.065" 






1.2 




HHIOO 


2.4 




3087 


0.032 (0.010) 


0.041 (0.005) 


1470 


0.15 (0.06) 


0.21 (0.09) 


8(3) 


NGC7538:IRS1 


3 




3065 


0.078 (0.013) 


0.052 (0.014) 


1455 


0.32"= 


0.24 (0.04) 


13 



° Water column density obtained from the integrated depth of the 3 /im band: NGC7538:IRS9/Allamandola et al. 1992; 

AFGL7009S/Dartois et al. 1999a (from 6 pra feature); W33A/Gibb et al. 2000 (obtained by fitting the 3 /im feature with 

laboratory analog spectra); AFGL 989/Smith et al. 1989; AFGL 2136/Schutte et al. 1996a; Ehas 29/Boogert et al. 2000; 

S140:IRSl(or GL2884)/Willner et al. 1982; W3:IRS5/Allamandola et al. 1992; MonR2:IRS3/Smith et al. 1989; HH100(or 

RCra:IRSl)/Whittet et al. 1996; NGC7538:IRS1/Willner et al. 1982. 
' Ratio of flux at 45 and 100 /im; from Keane et al. 2001 

Ground-based observations (Sellgren et al. 1995, Brooke et al. 1996, Brooke et al. 1999), unless otherwise noted. Error noted 

in parentheses. 

Keane et al. 2001; error typically 20 %, unless otherwise noted 
" From the ISO-SWS database 

^ It is unclear whether the 3 /im feature of MonR2:IRS3 is fully caused by H2O ice, due to the large width of this band (Smith 

et al. 1989). Therefore this number should be considered an upper limit. 
" Observed by ISO-SWS (Bregman et al. 2001) 

Boldface entries give the 3.26 /im depth derived from ISO-SWS data (Bregman et al. 2001). Errors are given in parentheses. 
' The uncertainty, due to the error in the bandstrength, is a factor ~ 2. 



al. 2000) extend into this region. Indeed, by excluding it 
we obtain a smoother baseline than Gibb et al. (i.e., third, 
rather than fourth order polynomial). 

Although there is a slight offset between the SWSl and 
SWS6 data, their overall spectral shape is quite consistent. 
Therefore the SWS6 data were converted to an optical 
depth scale by subtraction of the baseline of Fig. 0. The 
result is shown in Fig. |[ The small offset between the 
SWSl and SWS6 data was compensated by normalizing 
to in the 1900 - 1800 cm'^ region. 

The CH3 deformation feature of methanol falls near 
6.85 /im and will contribute a fraction of the interstel- 
lar band. The methanol column density towards W33A, 
as determined from the vj, and combination bands near 
3.54 and 3.9 /im, respectively, equals 1.85 x 10^^ cm~^ 
(Dartois et al. 1999b). Fig. || compares the optical depth 
spectrum of W33A with the band of methanol scaled 
to this column density. This feature was measured in a 
mixture H2O/CH3OH/CO2 = 1/1/1 (Ehrenfreund et al. 
1999). This matrix should be representative of the ices 
in which methanol in circumstellar objects is embedded 
(Gerakines et al. 1999; Dartois et al. 1999a). It is clear 
from Fig. that methanol only gives a small contribu- 



tion to the observed 6.85 /im band of ~ 20 % (see also 
Grim et al. 1991; Schutte et al. 1996). To enable a more 
suitable comparison with the laboratory NH^ spectra, we 
have subtracted this small contribution from the W33A 
spectrum. The resulting curve is shown in Fig. ^ Since 
the correction is small, using spectra of methanol in other 
matrices, such as pure methanol, or a water-dominated 
ice, gives a very similar result. 

Fig. ^compares the 6.85 /im band of W33A (after sub- 
traction of the methanol feature) with the band of NH^ 
as produced in one of our water-dominated samples after 
photolysis and warm-up (expt. 4). To take into account the 
variety of dust temperatures probed by the line of sight, 
the laboratory spectrum is an average of the 12, 120, and 
180 K (x2) spectra. We note that good overall matches 
can be obtained with all photolysed water-dominated ices 
(expts. 3, 4 and 5), since they produce little substructure 
in the 1400 - 1250 cm^^ (7.1 - 8.0 /tm) range. On the 
other hand, the photolysed ice mixtures with less water 
show strong peaks in this region due to NOj^ and NO^ 
(Sect. 3.1; Fig. H), which have no interstellar counterparts. 

Due to the superposition of ices at different temper- 
atures along the line of sight, components of different 
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Fig. 9. The optical depth spectrum of W33A after cor- 
rection for CH3OH (Fig. H), compared to: a. Photolysed 
H2O/CO2/NH3/O2 = 10/2/1.1/1 (expt. 4), average of 12 
K, 120 K, 180 K (x2) spectra. Curve a has been offset for 
clarity. 

volatility will have different average temperature. This is 
illustrated by the temperature averaged laboratory match 
of the W33A spectrum (Fig. H). The NH^ band arises pri- 
marily in the 120 and 180 K components, since at 12 K 
the intensity of this feature is small (Fig. ||) . On the other 
hand, the H2O ice features are produced in the 12 and 120 
K spectra, since at T > 170 K H2O sublimes. Thus, the 
average temperature of the solid H2O in the laboratory 
spectrum is 66 K, while the average temperature of NHJ 
is ~ 150 K. This clearly shows that temperatures derived 
by fitting the 6.85 /im band should not be taken as repre- 
sentative of the temperature of more volatile components 
located in cooler regions along the line of sight. 

4.1.2. MonR2:IRS3 

MonR2:IRS3 is an extreme object in a number of ways. 
First, the solid CO2 abundance is exceptionally low (~ 
1 %; Keane et al. 2001). Second, its 6.85 /im feature is 
the reddest of aU YSO's (Table |). Third, the strength 
of the 6.85 /im feature relative to the column density of 
H2O is very high (Table |3|). This all may indicate that 
the ice along the line of sight is exceptionally warm, as is 
also indicated by the high flux ratio F(45 /im)/F(100 /im) 
(Keane et al. 2001; Table |^). Indeed, the large relative 
intensity of the 6.85 /im band, assuming that the carrier 
of this band is a relatively refractory species, suggests that 
H2O ice has begun to evaporate along most of the line- 
of-sight. This is confirmed by the large column density 
of water vapor (Boonman et al. 2000). Ice sublimation 
should greatly enhance the visibility of features of more 
refractory components, making MonR2:IRS3 an excellent 



Fig. 10. The optical depth spectrum of MonR2:IRS3 
(From ISO/SWS; Keane et al. 2001) compared with some 
photolysed laboratory samples: a. II2O/CO2/NII3/O2 = 
1/0.5/0.7/1 at 160 K (expt. 2); b.same, at 170 K; c. 
10/2/1.1/1 at 180 K (expt. 4); d. 10/1.2/0.9/0.9 at 180 
K (expt. 5). The laboratory spectra are offset for clarity. 



testing ground for the presence of NII4 and associated 
counterions. 

Fig. |l^ compares the 6.85 /xm band of MonR2:IRS3 
with the NH|^ features obtained in several of our ex- 
periments. This set of experimental data will further down 
be used to evaluate the presence of other interstellar NHJ 
bands. No correction was made for the contribution of 
methanol to the 6.85 /im band, since no CII3OII has been 
observed towards this source (abundance < 4.5 %; Dartois 
et al. 1999b). The best matches are obtained for the pho- 
tolysed water-rich ices (expts. 4 and 5), since the water- 
poor ice (expt. 2) gives strong features at 1385 and 1335 
cm~^ due to NO3" which are not observed. In all cases, the 
position of the laboratory band is somewhat blue-shifted 
from the MonR2:IRS3 feature. This situation could be 
mended by selecting laboratory spectra at higher temper- 
atures. However, for such high temperatures a broad band 
near 1320 cm~^, which is caused by a superposition of fea- 
tures of various negative ions (Fig. ||, Sect. 3.1) shows up 
strongly, resulting in a decrease of the overall quality of 
the match. Even at the temperatures selected for Fig. [l^ 
too much absorption appears to be present in the 1350 
cm~^ region. However, this difference could be caused by 
the original choice of baseline, which made the a priori 
assumption of little absorption in the 1350 cm^^ region 
(Keane et al. 2001). As for W33A, considerable absorption 
may be present in this region. Thus, the position of the 
6.85 fim feature of MonR2:IRS3 is not fully reproduced by 
the experimental spectra, although the difference is small 
(~ 3-10 cm^^) for the three H2 0-dominated samples. In 
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Fig. 11. Position of the 6.85 iim feature vs. F(45 
Aim) /F( 100 urn) for the YSO's hsted in Table |. 



view of the strong matrix dependence of the position of 
this band (Fig. |3|), it seems plausible that the offset may 
be caused by a shift associated with the composition of 
the interstellar matrix. 

Besides the 6.85 /im band, the 6 /im feature of 
MonR2:IRS3 also has a counterpart in the laboratory 
spectra. Especially expt. 4 at 180 K gives a band which 
is quite similar (Fig. 



10) 



At this temperature most H2O 
has sublimed, and the laboratory features near 6 /im are 
dominated by a superposition of bands from various nega- 
tive ions, i.e., HCO3" and one or more unidentified species 
(Fig. ^; Sect. 3.1). This suggests that most of the 6 fiui 
feature of MonR2:IRS3 is caused by non-volatile compo- 
nents rather than H2O ice. Such a conclusion is consistent 
with the observations, which show that the 6 /xm feature 
of MonR2:IRS3 is 3.5 times stronger than expected at the 
H2O column density derived from the 3 /xm band (Keane 
et al. 2001). 



conclusion was already formulated by Keane et al., with 
the difference that they interpreted the change as due to a 
systematic variation with line of sight temperature in the 
contribution of two independent components of the 6.85 
fj,m feature. It must be stressed that interpretation of the 
temperature dependence of the band as either a variation 
of independent components, or a shift of a single compo- 
nent, is solely a matter of interpretation and cannot be 
decided on the basis of the interstellar band profile, since 
the feature does not show substructure. 

The systematic redshift of the 6.85 /im band with line- 
of-sight dust temperature is in good agreement with the 
trend observed for the 1^4 feature of the ammonium ion. 
Comparing Fig. ^ and Fig. ^ of Sect. 3.2.1, it can be 
seen that the spectral region in which the interstellar fea- 
ture is found falls well within the range spanned by the 
position of the band. Thus, the systematic redshift of 
the interstellar 6.85 //m band gives strong support to its 
identification with the 1^4 mode of NH^ . 

The range of positions of the interstellar band corre- 
sponds to laboratory temperatures of 120 - 240 K. It 
must be noted though that the position of the feature 
not only depends on temperature, but also on the ice ma- 
trix (Fig. Therefore, this temperature range should be 
taken as a rough indication. Nevertheless, it seems clear 
that solid NH^ is found at higher temperatures than more 
volatile ice components. Matching the profile of the 3 /im 
H2O feature typically yields 20 - 80 K (Smith et al. 1989), 
while matching the profile of the 15.2 /im bending mode 
of CO2 gives 110 - 140 K (Gerakines et al. 1999). This 
difference in temperature is consistent with the refractory 
nature of NH^ , which will survive in the warmer regions 
close to the embedded source where the more volatile ice 
components have sublimed. In addition, the band strength 
of NH^ increases considerably above ~ 100 K (Sect. 3.1; 
Fig. H). This effect will minimize the contribution of the 
coldest regions to the 6.85 /im band. 



4.2. The redshift with temperature 

It was shown in Sect. 3.1 that the ammonium band shifts 
strongly towards the red with increasing temperature. 
Even though the position of the NH^ band depends con- 
siderably on the starting matrix, this general trend held 
for all our experiments. Thus, if the ammonium ion is the 
carrier of the interstellar 6.85 /im feature, a correlation 
should be found between the position of the 6.85 /im fea- 
ture and indicators of the dust temperature along the line 
of sight. 

Fig. In] plots the position of the 6.85 /im band in the 
data set of Keane et al. (2001) against the ratio of the 
fluxes at 45 and 100 /xm as measured by ISO/LWS. This 
ratio is a general indicator of the dust temperature along 
the line of sight. It can be seen that there is a clear corre- 
lation between the position of the band and the flux ratio, 
in the sense that the interstellar band shows a strong red- 
shift with increasing average dust temperature. A similar 



4.3. Additional ISfBj features 

Besides the feature near 6.85 /im, the ammonium ion 
shows a number of additional features, i.e., 2i/2, 1^1 +1^5-, 
and 21^4. These are centered near ~ 3200 cm~^ (3.12 /im), 
~ 3060 cm-i (3.27 /im), and ~ 2860 cm"! (3.50 /im), 
respectively (Figs. |^ & H, Sect. 3.2.1). The optical depth 
of each of these bands is typically 30 % of that of the 1^4 
feature. 

The 2f2 feature would be very difhcult to distinguish, 
since it nearly coincides with the peak of the 3 /im water 
band (Smith et al. 1989). The two other bands, while rela- 
tively weak, may be observable as substructure on the red 
wing of the water band. Indeed, features have been ob- 
served towards many protostellar sources near these posi- 
tions (e.g., Sellgren et al. 1995; Brooke et al. 1996; Brooke 
et al. 1999; summary in Table |). For example. Fig. [l| 
shows the mid-IR spectrum of MonR2:IRS3 obtained by 
ISO-SWS (Gibb et al. 2001) and by ground-based spec- 
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Fig. 12. Comparison between the 3.26, 3.48 and 6.85 [ivn features of MonR2:IRS3 and the V2 + i^4, 2z^4 and features 
of NH^ in a number of laboratory samples. The various NH^ features are indicated by vertical dashed lines. The 
labelling of the laboratory spectra corresponds to the description given with Fig. Dotted lines indicate the adopted 
baselines for the two features. The observational data comprise ground-based observations between 3.16 - 3.65 /xm 
(Sellgren et al. 1995; bold solid line) and ISO/SWS data (2.6 - 8 ^im; Gibb et al. 2001; soHd Hue). The MonR2:IRS3 
data were converted to optical depth by subtraction of a 2"'' order polynomial baseline defined by fitting the continuum 
regions around 3800 and 2000 cm^^ in the original flux data. All spectra were smoothed to resolution 400. The spectra 
have been offset for clarity. 



troscopy (between 3.16 - 3.65 ^m; Sellgren et al. 1995). 
As already mentioned above, this source appears to be 
exceptionally warm and therefore should have enhanced 
abundances of refractory components like NH4" . 

To facilitate comparison with the laboratory data, a 
continuum has been subtracted from the MonR2:IRS3 
original flux data (Ground-based as well as ISO-SWS). 
This continuum was obtained using the ISO-SWS data 
by drawing a second order polynomial through the zero 
absorption regions around 3800 and 2000 cm~^. Before 
continuum subtraction, the ground-based data were mul- 
tiplied by 0.88 to match them to the ISO-SWS spectrum. 

A small excess is apparent in the ISO-SWS data rela- 
tive to the ground based spectrum around 3.29 /im. This 
excess is obviously caused by the presence of the well- 
known emission feature in the wide aperture ISO data 



(This is a common phenomenon for ISO observations to- 
wards YSO's; Gibb et al. 2001). 

As already discussed by Sellgren et al., MonR2:IRS3 
shows two features at 3.26 and 3.48 ^m. Fig. |lj also shows 
the V2 -I- 1^4 and Iv/^^ NH^ bands as obtained in our ex- 
periments. To allow a detailed comparison, we subtracted 
baselines from both the laboratory and interstellar spec- 
tra, as indicated in Fig. 02. The resulting optical depth 
curves are shown in Fig. Il3. Only the ground-based ob- 
servations were used, due to the "contamination" of the 
ISO spectra with the 3.3 /zm emission feature. 

The comparison in Figs. ^ and |l^ indicates a good 
general agreement between the 3.26 and 3.48 /im features 
of MonR2:IRS3 and the V2 + and 2v^ NH| bands. 
Nevertheless, no single spectrum is able to fit in detail 
both the position and relative intensity of the observed 
features. However, the variation among the laboratory 



12 



W. A. Schutte et al.: NH4 and the 6.85 fim band near Young Stellar Objects 



3.3 3.4 



3.5 3.6 



-0.12 



g- -0.04 
Q 



D- 0.00 



- 


1 ' 1 ' 1 ' 


1 ' 1 ' 1 
















\ ^"^^ //\ 






\ / r-^ V ^ 






1 ' 1 ' 1 


/MonR2:IRS3 
' 1 ' 



3200 



3100 



3000 2900 

V (cm" ) 



2800 



2700 



Fig. 13. Comparison between the 3.26 and 3.48 /im (3075 
and 2875 cm^^) features of MonR2:IRS3 and the P2 + 
and 21/4 bands of NHJ in the laboratory samples of Fig. 
after subtraction of the baselines indicated in this figure. 
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spectra shows that the features are quite matrix and tem- 
perature dependent, and it therefore seems plausible that 
this moderate difference could derive from matrix effects. 
Besides the agreement in profile, the agreement between 
the intensity of the laboratory and the MonR2:fRS3 fea- 
tures lends support to an identification with NH^. 

It is clear from Fig. |lj that if the 3.26 and 3.48 /im 
features of MonR2:IRS3 are caused by the ion would 

also partially account for the overall broad absorption in 
the 3.1 /im feature. Indeed, some 50 % of the unidentified 
excess absorption between 3.25 and 3.8 /im (3080 - 2630 
/im~^), which cannot be accounted for by H2O ice (Smith 
et al. 1989), could arise from the ion. 

If the 3.26, 3.48, and 6.85 /im features are all due to 
the same carrier, a tight correlation between these bands 
is expected. Ground based observations of the 3.26 /im 
feature exist for four of the objects whose 6.85 /im feature 
has been observed by the SWS. (Table ||). For these ob- 
jects the data give a stable ratio t(3.26)/t(6.85) — 0.22 ± 
0.02. The constancy of the ratio is even somewhat better 
than what could be expected on the basis of the rather 
large uncertainty in the depth of the 3.26 /im band (Table 
4). Besides from telluric interference, this error derives 
from the baseline definition, since this shallow feature is 
located in a complex spectral region, where various spec- 
tral features are located, such as the broad feature of am- 
monia hydrate around 3.5 /xm (2860 cm^^) and the 3.36 
/im (2980 cm~^) feature of CII3OH, as well as possibly 
structure caused by scattering (Dartois & d'Hendecourt 
2001; Brooke et al. 1999). These features may contribute 
to the broad long wavelength shoulder of the 3 /im band, 
on which the feature is superimposed. 



For some of these objects ISO-SWS observations on 
the 3.26 /im feature have also been reported (Bregman et 
al. 2001). Using these data, there is little or no correlation 
between the 3.26 and 6.85 /im bands (Table |). Indeed, 
the ground based and ISO data give conflicting results on 
the depth of the 3.26 /im band, most notably in the case 
of p Oph Elias 29, where the discrepancy is at least a fac- 
tor of 3. It not quite clear what causes these differences 
and which dataset is most reliable. For Elias 29, Bregman 
et al. used a speed 3 scan in AOT SWSl. However, a 
higher quality SWS6 scan in this region docs not show 
the 3.26 /im feature, consistent with the ground-based ob- 
servations (Boogert et al. 2000; Boogert, private commu- 
nication). A general problem with the ISO data could be 
their wide aperture. For this reason, the 3.3 /im emission 
feature often shows up in the SWS spectra of embedded 
YSO's (Gibb et al. 2001; see also Fig. Indeed, some 
3.3 /im emission seems consistent with the ISO spectra of 
S140 and GL 989 (Gibb et al. 2001). The presence of even 
a small amount of this feature would severely complicate a 
correct baseline definition and therefore introduce a large 
error in the depth of the shallow 3.26 /im band. On the 
other hand, in ground-based observations the center of the 
3.26 /im band has an enhanced noise level due to telluric 
methane absorption. However, the ISO and ground-based 
observations also deviate for the 3.48 /im band (Table 4), 
even though it falls in a relatively undisturbed region. This 
deviation therefore reinforces the concern that errors in 
the baseline definition caused by the 3.3 /im emission fea- 
ture is a major pitfall when using ISO-SWS data to study 
the 3.26 and 3.48 /im absorption bands. 

Good candidates to search the 3.26 /im feature are 
NGC7538:IRS9 and W3:IRS5, for which we predict, from 
the strength of the 6.85 /im band, t(3.26) w 0.064 and 
0.055, respectively. A detection of these bands would be 
an important test of the NH^ assignment. 

No tight correlation exists between the 3.26 /im and 
3.48 /im features. Relative intensities vary between 0.13 
(MonR2:IRS3) and 1.3 (Ehas 29; see Table |). While the 
relative intensity of the feature of MonR2:IRS3 is consis- 
tent with an assignment of the 3.48 /im band to NH^ 
(Fig. 13), it is clear that in general other species must 
also contribute to this feature. This is not too surprising, 
since additional structure is present in this region. In par- 
ticular, the long wavelength shoulder or wing of the 3 /im 
H2O band, which extends from - 3.1 - 3.8 /im (3200 - 2600 
cm^^; Smith et al. 1989; their Fig. 6) produces structure 
in this region (Smith et al. 1989; their Fig. 4). This ab- 
sorption will blend in with any additional absorbers, e.g., 
the CH stretching mode of aliphatic hydrocarbons and 
the ^ 3.5 /im feature of ammonia hydrate (e.g., Dartois 
& d'Hendecourt 2001). 



4.4. The 3.26 jim feature and PAHs 

Previously, an assignment of the 3.26 /im feature to the 
CH stretching modes of polycyclic aromatic hydrocarbons 
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has been proposed (PAHs; Sellgren et al. 1995; Bregman 
et al. 2000). Such an assignment would of course Umit 
the contribution of the 1/2 + 1^4 NHj feature to this band 
(Sect. 4.3). Here we quantitatively assess the contribution 
of the PAH CH stretching mode to the 3.26 /im feature of 
MonR2:IRS3. 

In support of, at least, part of the 3.26 /^m feature 
deriving from the PAH CH tretching mode is the obser- 
vation of a weak absorption band at 11.2 /im towards 
MonR2:IRS3, which is ascribed to the out-of-plane bend- 
ing mode of isolated CH groups on the aromatic skele- 
ton (Bregman et al. 2000). The intensity of this band can 
be used to constrain the PAH contribution to the 3.26 
/im feature. The integrated optical depth are ri„t(11.2) = 
0.81 cm~^ and ri„t (3.26) = 5.6 cm~^. The ratio of the 
band strength of the CH stretching mode to that of the 
CH bending mode for isolated CH groups equals 0.5 for 
gaseous, neutral PAHs. For ionized or solid state PAHs the 
ratio is at least 4 times smaller (Bregman et al. 2000 and 
references therein). Therefore the CH groups associated 
with the 11.2 /im feature can, at most, cause an absorp- 
tion of Ti„((3.26 /im;PAH) = 0.4 cm~^, i.e., only 7 % of the 
observed intensity. Besides isolated CH groups, aromatic 
species will have CH groups on their periphery which have 
1 or more adjacent CH groups. For these the CH bending 
mode falls at longer wavelength (11.6 - 13.6 /im; Hony et 
al. 2001). Besides discrete features near 12.7 and 13.5 /im, 
the bending modes of these groups give rise to a broad 
plateau between 11-13 /im caused by overlapping ab- 
sorptions. In emission, the fraction of the total emission 
associated with the CH bending modes (i.e., 11.2, 12.7, 
13.5 /im bands and the plateau) which is emitted by the 
11.2 /im band varies between 0.3 - 0.7 (Hony et al. 2001). 
Since the band strength of the CH bending mode for ad- 
jacent aromatic CH groups is two times weaker than for 
isolated CH groups, this implies that these groups could 
add Tint = 0.4 - 1.9 cm"-'^ to the 3.26 /im feature. This 
implies that 14 - 40 % of the 3.26 /im band may be due to 
the PAH CH stretch. This fraction could actually be much 
lower if a large fraction of the PAHs would be frozen on 
the grains, causing the band strength of the CH stretching 
mode to drop by a factor ~ 4 (Joblin et al. 1994). Indeed, 
it seems probable that in a high density environment like 
the circumstellar regions of a young stellar object a large 
fraction of highly refractory molecules such as PAHs would 
accrete on the dust grains. In conclusion, it seems prob- 
able that the contribution of PAHs to the 3.26 /im band 
is rather small. This agrees with our result in Sect. 4.3, 
namely, that, based on the assignment of the 6.85 fim band 
to NH|, most of the 3.26 /im feature should be caused by 
the NH^" 1^2 + va band. 

Contrary to our result it was concluded by Bregman et 
al. (2000) that PAHs may account for the entire 3.26 /im 
absorption band. Some important differences between his 
calculations and ours must be pointed out. First, the in- 
tegrated intensity of the 3.26 /im band was approximated 
by depth x FWHM in this earlier paper, while we use ex- 
act numerical integration, arriving at a 1.7 times higher 



intensity. Second, Bregman et al. estimated that typically 
only 17 - 30 % of the entire emission in the CH bend- 
ing region originates from the 11.2 /im band. However, 
a thorough recent survey of this region in a wide variety 
of objects (Hony et al. 2001) showed that this fraction is 
considerably higher, i.e., 30 - 70 %. 

4.5. The abundance o/"NH^ and the counterions 

Applying the intrinsic strength of the 1^4 NHj band, the 
NH^ column density can be derived from the observa- 
tions. With A(NH+) = 4.4 X 10"" cm ion^i (Sect. 3.3; 
uncertainty factor 2), we find for W33A, with an inte- 
grated optical depth of Ti„f(6.85 /im) = 57 cm~^ (after 
correction for the CH3OH contribution), N(NH^) = 1.3 
X 10^^ cm-2. For W33A we adopt N(H20) = 1.2 x 10^^ 
cm~^. This number was derived by Gibb et al. 2000 by 
fitting the blue wing of the 3 /xm band with laboratory 
ice bands. It must be noted that those authors adopted 
a slightly lower column density as their end result, based, 
in part, on matches of the blue wing of the 3 /im band of 
W33A with other interstellar 3 /im features. However, the 
correspondence of the W33A 3 /im band with the other 
interstellar features was considerably less than the corre- 
spondence with the selected laboratory ice bands (presum- 
ably caused by source to source variations in ice composi- 
tion and temperature). Therefore, we prefer to adopt here 
the value which is solely based on the laboratory compar- 
ison. This gives an NH| abundance of 11 % for W33A, 
with an uncertainty of a factor 2 due to the bandstrength 
error. 

From the data summarized in Table ^, it can be de- 
rived that the NH^ abundance varies from 3 % to 17 % 
relative to H2O. Here we use a constant correction for the 
contribution of methanol to the 6.85 /im band of 20 %, as 
for W33A. This is perhaps a bit large, since in general the 
methanol abundance in YSO's is lower than for W33A 
(Dartois et al. 1999b). However, the correction is small 
anyway, and no attempt was made for a detailed source 
to source analysis. One exception is GL7009S, which has 
an exceptionally high abundance of solid methanol, for 
which we adopted a 30 % correction of the 6.85 fim fea- 
ture (Dartois et al. 1999a). 

Clearly, an equal amount of negative charge should be 
present in the ices to balance NH^. Recently, Gibb et al. 
(2000) investigated the full inventory of ice species along 
the line-of-sight towards the high mass YSO W33A. Their 
analysis showed that the negative species that have been 
detected towards this source, i.e., OCN~ and HCOO~ (as- 
suming that OCN~ is the carrier of the interstellar XCN 
feature; Grim & Greenberg 1987; Schutte & Greenberg 
1997; Demyk et al. 1998; Hudson & Moore 2000; 2001; 
Novozamsky et al. 2001) have abundances of 3.2 % and 
0.8 %, respectively (these values are slightly lower than 
those of Gibb et al., due to the 10 % higher H2O col- 
umn we adopt for W33A). Since the NH^ abundance for 
this source is ~11 %, this implies that only ~ 30 % of 
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the countercharge is provided by the observed negative 
species. Our experiments suggest that the residual charge 
could be provided by a combination of ions such as HCOj^ , 
NO^ , NO^ . As was shown, a mixture of these species em- 
bedded in an H20-dominated ice does not produce signif- 
icant spectral structure. Their total abundance would be 
~ 7 % of H2O for W33A. For other sources, the fraction of 
the charge balance that can be provided by the observed 
negative species (i.e., OCN~) is even smaller (Gibb et al. 
2000). Thus, assuming that they make up for the residual 
balance, the total abundance of "invisible" anions in the 
ices near YSO's would be in general - 70 - 100 % of NH|, 
i.e., 3 - 15 %. 

It has been argued that the poor correlation between 
the 4.62 /im feature of OCN~ and the 6.85 /im band dis- 
agrees with an identification of this feature with NH^ 
(Keane et al. 2001; cf. Figs. & |l^. However, as pointed 
out above, OCN~ provides only a small fraction (< 20 %) 
of the countercharge. In such clear correlation be- 

tween the OCN~ and NH4 abundances is only expected if 
the make-up of the anion mixture is stable between YSO's. 
The lack of correlation shows that the fraction of OCN~ 
in the total ensemble of anions strongly varies. 

5. Astrophysical Implications 

When the NH^ assignment was originally made (Grim et 
al. 1989a), the identity of the counterions was uncertain. 
This was widely considered a prime concern (Schutte et al. 
1996a; Keane et al. 2001; Demyk et al. 1998; Tielcns and 
Whittet 1997). Our new experiments demonstrate that in 
II2O dominated ices the charge balance can be achieved 
without introducing spectral structure of negative ions 
that is inconsistent with the observations. Thus it appears 
that the objections to the NH^' assignment have been met. 
All things considered, NH^ now stands as the best assign- 
ment for the interstellar 6.85 /im feature (apart form the 
generally small fraction caused by CH3OH; Sect. 4), based 
on the band's position, shape, and width. Moreover, the 
ammonium ion matches both the 6.85 and 3.26 /im fea- 
tures, and is able to reproduce the temperature-dependent 
redshift of the 6.85 /im band. Chemical considerations also 
support the assignment as NHj is easily produced by ei- 
ther photolysis or ion irradiation of known, or suspected, 
interstellar molecules, or by mild warm-up of ammonia 
and acids which could be present in the ice (see below). 
A straightforward observational test would be to probe 
the tight correlation between the 3.26 and the 6.85 /im 
features which is expected if both are due to NII4" . 

It is unclear whether interstellar NH^ could arise from 
irradiation or photolysis as in our experiments. The weak- 
ening of the 3.4 fim CH stretching mode of aliphatic grain 
material in goig from diffuse to dense regions indicates 
that energetic processing could play an important role 
in dense regions (Mufioz Caro et al. 2001; Mennella et 
al. 2001). However, the processing of our samples gives 
rise to some weak features which have not yet been ob- 
served, most notably the O3 band at 1038 cm~^, and the 



unidentified 1270 cm~^ feature (Fig. ^). Neither of these 
bands are visible in the ISO-SWS observations of YSO's 
(Fig. Q Gibb et al. 2001). Alternatively, the formation 
of NHJ may not be associated with energetic process- 
ing at all. Perhaps NH3 and acids like H2CO3, HNO3, 
HNO2, HNCO or even H2SO3 and H2SO4 can be formed 
by surface chemistry. Whether such relatively large species 
are formed by surface reactions depends on the extent at 
which the exothermicity of the reactions can be channeled 
into translational energy of the products (Brown 1990; 
Hasegawa et al. 1992). Acid- base reaction and ion forma- 
tion could also be part of the surface reaction sequence, or 
be promoted by warm-up or a small amount of radiation 
(Novozamsky et al. 2001, van Broekhuizen et al. 2003). 

While in this view the anions would be undetectable 
in the solid state, possibly the corresponding acids could 
be found in the gas phase in the hot cores of star form- 
ing regions where the ice mantles evaporate. A search for 
molecules like H2CO3, HNO3, HNO2 or H2SO3 in such 
regions seems therefore warranted. 

From the limited sample of objects listed in Table || it is 
unclear whether NH4' has different abundances for low and 
high-mass objects. The low-mass protostar Elias 29 in the 
p Oph cloud has an NH^ abundance of 3 %, substantially 
less than for the high mass sources. However, the other 
low mass object in our sample, HHIOO, gives a relatively 
high abundance of 8 % (albeit with a large errorbar). More 
observations of the 6.85 /im feature of low mass protostars 
and background field stars are desirable to clarify which 
dense environments favor the formation of NH^ . 

In hot core regions where elevated grain temperatures 
cause sublimation of the icy mantles, large gas phase abun- 
dances of ammonia have been found ((1 - 10) x 10~^ rel- 
ative to hydrogen; Blake et al. 1987; Heaton et al. 1989; 
Cesaroni et al. 1994). This considerably exceeds the abun- 
dance of NH3 in cold dense regions (N(NH3)/N(H2) « 
10"^; Federman et al. 1990), indicating that ice sublima- 
tion is an important source of NH3. However, recent obser- 
vations failed to detect the 2.21 /tm feature of NH3 towards 
W33A, indicating that the NH3 abundance in interstellar 
ice may be low (< 5 %; Taban et al. 2003). A high abun- 
dance of solid NH^ in interstellar ices could alleviate this 
dilemma, since sublimation of this compound takes place 
during warm-up upon a reverse acid-base reaction and the 
re- formation of NH3. 

NH^ can only be formed from NH3. The presence of 
solid NH^ thus indicates that solid NH3 was efficiently 
formed during the ice accretion, presumably by hydro- 
genation of atomic nitrogen on the grain surface (Hiraoka 
et al. 1995). Whether in dense regions gaseous nitrogen 
will reside in its molecular or its atomic form depends on 
the depletion of oxygen, since reactions of atomic nitro- 
gen with OH is an essential step towards N2 formation 
(Charnley & Rodgers 2002). Thus the quantity of NH| 
that is present in the ices gives important information on 
the gas phase conditions at the epoch of condensation. 
Thus probing the strength of the 6.85 /im band could be- 
come an invaluable tool for studying the chemical condi- 
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tions associated with the different regions in dense clouds, 
i.e., quiescent, as well as near embedded high mass and 
low mass YSO's. 



6. Conclusions 

Experiments involving UV photolysis of 
H2O/CO2/NH3/O2 ice mixtures efficiently produce 
ammonium (NH4 ) together with probably HCO3 , NO^, 
and NO2 • These ions are produced by acid-base reactions 
between the NH3 and photochemically produced acids. 
The Pi mode of NH4 gives a strong feature at 6.85 
jjLva.. Other NH4 features fall at 3.27 and 3.50 ^m. Due 
to strong interaction of the ion with the matrix, the 

NH4 band shifts strongly redward with increasing 
temperature. The infrared signature of the counterions is 
very weak for H20-dominated ices. 

Comparison to the 6.85 /im feature towards embed- 
ded young stellar objects shows that the 1^4 NH4 feature 
obtained in our experiments provides a good match. The 
observed red shift of the interstellar band with increasing 
line-of-sight dust temperature is also consistent with this 
assignment. Furthermore, features are observed at 3.26 
and 3.48 /im towards these objects which could be (partly) 
due to the ammonium ion. The implied abundance of NH4" 
relative to H2O is 3 - 17 %. 

The two anions observed in interstellar ices, OCN~ 
and (possibly) HCOO~ , are insufficiently abundant to bal- 
ance the positive charge from the NH4' present. Therefore 
heavier negative ions of the kind formed in our experi- 
ments appear to be required. It is yet unclear whether 
such ions are formed by energetic processing, as in the 
laboratory experiments, or by grain surface chemistry. 

Future laboratory studies should include analysis by 
mass spectroscopy to exactly define the nature of the pho- 
tochemically produced acids. Sulfur could be added to the 
mixture to extend the number of acids that are produced. 
On the observational side, the NH^ assignment could be 
tested by probing the correlation between the 3.26 and 
6.85 /J,m bands. In addition, it woidd be very interesting to 
search for rotational emission bands of evaporating acids 
like HNO2, H2CO3 or H2SO3 in hot core regions. 
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